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ABSTRACT 

NASA’s re-purposed Kepler mission - dubbed K2 - has brought new scientific opportunities that were not 
anticipated for the original Kepler mission. One science goal that makes optimal use of K2’s capabilities, in 
particular its 360-degree ecliptic field of view, is galactic archaeology - the study of the evolution of the Galaxy 
from the fossil stellar record. The thrust of this research is to exploit high-precision, time-resolved photometry 
from K2 in order to detect oscillations in red giant stars. This asteroseismic information can provide estimates 
of stellar radius (hence distance), mass and age of vast numbers of stars across the Galaxy. Here we present the 
initial analysis of a subset of red giants, observed towards the North Galactic Gap, during the mission’s first 
full science campaign. We investigate the feasibility of using K2 data for detecting oscillations in red giants 
that span a range in apparent magnitude and evolutionary state (hence intrinsic luminosity). We demonstrate 
that oscillations are detectable for essentially all cool giants within the logp range ^ 1.9-3.2. Our detection is 
complete down to Kp ^ 14.5, which results in a seismic sample with little or no detection bias. This sample 
is ideally suited to stellar population studies that seek to investigate potential shortcomings of contemporary 
Galaxy models. 

Subject headings: stars: fundamental parameters — stars: oscillations — stars: interiors 


1. INTRODUCTION 

The study of red giant stars has arguably been one 
of the greatest success stor ies of NASA’s Kepler mission 
(e.g., iGarcfa & Stelloll2015L and references herein). How¬ 
ever, a failure of the second of four momentum wheels ended 
the mission in 2013 because the spacecraft could no longer 
acquire stable pointing towards its original field of view. For¬ 
tunately, ingenious use of the remaining spacecraft capabili- 
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ties by NASA and Ball Aerospace engineers rejuvenated the 
mission as K2 - a mission capable of stable pointing at any 
field along th e ecliptic for up to approximately three months 
per pointing (iHowell et al.ll201^ . In this configuration, the 
Kepler roll angle drifts due to a torque applied by solar radi¬ 
ation pressure, but this can be counteracted by thruster firings 
every six hours to maintain the spacecraft pointing. The K2 
mission has en abled a broad range of new science including 
stellar clus ters (INardiello et al.ll2 Q15h. nlanets around bright 
cool stars (ICrossf ield et al.ll2015l: I Sanchis-Oieda et al]|2015l: 


IVanderb urg et al.l 20151: Montet et alJl2015h, solar system ob¬ 

jects d^abo et al.l l2015h . s tellar activity (iRamsav & Dovlel 
(2015h . e clipsing binaries ( Conrov et alJ l2014h. a s teroseis- 
mology (iJefferv & Ramsa3l 20141 : iLund et al.ll2015l : IChaplinI 
IMll and, in particular, asteroseismological studies of the 
Galaxy. 

The potential for asteroseismic investigations of large pop¬ 
ulations of red giants aimed at Galactic studies was re- 
cently demonstrated using data from CoRoT and Kepler 
(iMiglio et al.l 120091 : IChaplin et al.l 1201 ll : iMiglio et'aP 120131 : 
ICasagrandell20 1 5h However, the scope of these early studies 
was limited for two reasons: the small number of distinct di¬ 
rection in the Galaxy probed by those missions, and the highly 
complex (and, at some level, not fully documented) selection 
function of the observed red giants (Sharma et al. in prepa¬ 
ration). With K2’s 360-degree coverage of the ecliptic the 
collated efforts from the K2 observing campaigns, provide a 
unique opportunity to probe different regions of the Galaxy, 
including the thin and thick disks, the halo, and the bulge, 
based on a purpose-built selection approach suitable for pop¬ 
ulation studies. 

In this Letter, we present initial results from the K2 Cam¬ 
paign 1 data. Based on a sample of red giants specifically 
selected to study stellar populations on a galactic scale, our 
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aim is to determine if we can obtain an unbiased sample of 
stars showing oscillations - a crucial first step if these stars 
are to be used for galactic archaeology studies. 


2. OBSERVATIONS AND LIGHT CURVE PREPARATION 

We used observations obtained as part of the K2 Galactic 
Archaeology Program Campaign 1 (Cl proposal 001059)'^. 
For the purpose of this Letter, we focused on a targeted set of 
red giant candidates based on their sp ectroscopic log q value s 
(logp < 3.8) drawn from APOGEE (iMaiewsk i et ai.ll2 010t). 
part o f the Sloan Digital Sky Survey 111 (lEisenstein et all 
1201 Ih . The APOGEE survey used a wide -held multi-object 
H-band spectrograph (I Wilson et al.l 1201011 on the 2.5-meter 
Sloan Eoundation telescope at Apache Point Observatory 
dGunn et alJl2006l) . The red gian ts discussed here wer e ob¬ 
served as part of the main survey (IZasowski et alJl2()T^ . The 
spectra were reduced, wavelength -calibrated, and co-added as 
described bv iNidever et al.l (l2015h . The determination of stel¬ 
lar parameters from the the automated pipeline is described 
by Garcia-Perez et al, in p reparation). We u sed the parame¬ 
ters from Data Release 12 (lAlam et al.ll2015l) : the calibration 
and verification of the APOGEE DR 12 results is described by 
iHoltzman et alJ (l2015h . Of the 121 red giant candidates that 
satisfy the above-mentioned logp cut, 117 were observed in 
the K2 campaign. The stars span 2.8 dex in logp and seven 
magnitudes in apparent magnitude (9 < Kp < 16), and hence 
serve as a suitable benchmark set to characterize the K2 data 
fidelity, including the ability to detect oscillations for different 
levels of intrinsic and apparent brightness. 

The photometric time series (light curves) from the raw 
K2 pixel data are currently n ot available from NASA . How¬ 
ever, light curves created by IVanderburg & JohnsonI (12014l) 
(VJ hereafter) are publicly available'* and we used these data 
in our study. Power spectra of the light curves denoted ‘cor¬ 
rected’ by VJ are shown in Fig. [T] (black curves) for three 
stars. The 80-day duration of the light curves provides a 
nominal frequency resolution of ^ 0.14 p,Hz, which imposes 
a lower limit on the frequency separation between overtone 
modes, Av, that we can reliably determine dHuber et al.l2010l : 
iHekker et alJl2012)). Through th e tight correlation between 
Av and i^max (IStello et al.ll2009li . this translates into a lower 
limit of about 10p,Hz in Vniax, and hence to the frequency 
range in which we can fully characterise the oscillations. By 
restricting our attention to higher-frequency stars, we are then 
free to apply a high-pass filter without compromising the os¬ 
cillation signal. We chose a boxcar filter with a width of two 
days, resulting in a cut-off frequency of about 3pHz (Fig.[Tl 
red spectra), without affecting the oscillation signal of our tar¬ 
get stars. Applying a high-pass filter gives a significant re¬ 
duction in the noise floor at all frequencies because it reduces 
spectral leakage of low-frequency power to higher frequen¬ 
cies, as seen by comparing the black and the red spectra in 

Fig.n 

Due to the roughly 6-hour drift and repointing cycle, K2 
data generally show trends on that time scale, and typically 
one d ata point is flagged unsuitab le during the repointing 
(see IVanderburg & JohnsonI l20f^ . The resulting regular 
gaps, combined with the slow drifts, can result in significant 
leakage of power towards higher frequencies in the power 
spectra. This offset can be avoide d if the gaps are filled 
(iGarcfa et al.ll2014l : iPires et aTll2015h . for which we used lin- 
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Eig. 1.— Power spectra of three suspected red giants. Their EPIC IDs 
(Huber et al. 2015, in preparation) are shown in each panel. An'ows mark the 
location of excess power assumed to be from stellar oscillations. Dotted lines 
are located at multiples of 47.2281^Hz, corresponding to the spacecraft roll 
angle re-positioning frequency. Spectra in black are based on Vi’s ‘corrected’ 
light curves. Red spectra show the result after high-pass filtering the light 
curves. Spectra in green include both a high-pass filter and an interpolation to 
fill short regular gaps in the data caused by the repositioning of the spacecraft. 


ear interpolation for all gaps of up to three consecutive data 
points. Larger gaps were not filled. The interpolation used 
only the two points tha t bracketed each g ap. In agreement 
with lGarcla et al.l(^2014^ : |Pires et al.l (120151) . we found that the 
stellar signal was not markedly affected. However, our ability 
to detect the oscillations was increased, particularly for stars 
oscillating at relatively high frequencies. This improvement 
is illustrated by comparing the red and green spectra in Fig.[T] 
While gap filling should be used with caution, it is generally 
safe when the gaps are short and few co mpared to the total 
number of data p oints in the light curve (iGarcla et al.ll2014t 
iPires et aTl 120151) . Here, we filled 3^% of all data points. 
For these K2 data, the gap hlling results in almost complete 
removal of peaks in the power spectra at the re-positioning 
frequency and its harmonics (Fig.[T] dotted lines), which are 
otherwise disruptive for the automated detection of oscilla¬ 
tions and extraction of the global seismic properties. In the 
following analysis, we use the high-pass hltered, gap-filled 
light curves. 

In future work, we will explore other scheme s for directly 
model ling instrumental variability. In particular, lAngus et al.l 
(1201 5l) have recently proposed a scheme that alleviates the 
need for gap filling in time series data to produce K2 power 
spectra that are less sensitive to systematic effects. 

3. OSCILLATION ANALYSIS 


” http://keplerscience.arc.nasa.gov/K2/index.shtml The power spectra of many Stars in our sample reveal clear 

https://www.cfa.harvard.eduravanderb/k2.html oscillations ranging from low-luminosity giants near the bot- 
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Fig. 2.— Power spectra of selected red giants representing the range in 
i^max, detectable from K2 long-cadence data. Their EPIC ID’s are indicated. 
The stars are ordered from most luminous (lowest oscillation frequencies of 
about 10/iHz) at the top, to the least luminous (highest oscillation frequency 
of about 270/iHz) at the bottom. The frequency range dominated by the 
oscillations is shown in black. 


tom of the red giant branch to stars more luminous than the 
red clump. The detected oscillations cover the same fre- 
quency range as e arly Kepler light curves of similar length 
(iHuber et al.ll2010|) . Fig.|2]presents power spectra of a repre¬ 
sentative subset of our sample. In each star, we see the os¬ 
cillation power excess forming a near regular series of peaks 
from overtone modes within a broad envelope, from which 
Vraax and /S.V can be measured. 

We carried out a syst ematic search for os cillations using the 
pipeline developed by iHuber et~aD (l2009ll . After inspecting 
the power spectra and the diagnostic output from the auto¬ 
mated detection algorithm, we classified the stars into three 
groups. A total of 59 stars provided clear detections of both 
Vraax and /S.V, 14 stars were classified marginal, and 44 were 


non-detections. Marginal detections refer to stars where ei¬ 
ther i.'niax or Av were not determined unambiguously. Non¬ 
detections are mostly those stars where we did not find any 
evidence of oscillations. Some were slowly pulsating (very 
luminous) giants, which in some cases did show evidence of 
oscillation power, but due to the small i^max and Av we are 
not confident in claiming those as marginal detections. Based 
on previous experience with Kepler data, most of the marginal 
detections and almost all of the non-detections extend over re¬ 
gions in log g-Tes space that render them unlikely to result in 
measurable and Ai/ values. 

To illustrate this property, we show logp and Teff from 
APOGEE in Eig. of all 117 stars in our sample superim¬ 
posed on a stellar evolutionary track (gray curve) of a 1.2Mfn, 
rough ly solar met allicity MESA mod el (iPaxton et alJ 1201U 
12013h taken from IStello et aP (l2013h . Prom stellar eyolu- 
tion models, we would generally not expect to find giants 
hotter than Teff = 5000 K (except for rare extremely metal- 
poor and/or massiye stars. Indeed, preyious results from long¬ 
term obseryations by Kepler, showed only two oscillating gi- 
ants hotter than APOGEE Tea = 5000 K (iPinsonneault et all 
12014 : one on the lower red giant branch and one red clump 
star out of 1916 stars in total. Our K2 results appear to be in 
line with those results, and we should therefore discount stars 
that appear to be hotter than this Tea threshold when assess¬ 
ing our detectability capabilities (Pig.[^, yertical dotted line). 
Howeyer, we do not adyice to apply a Tes selection for future 
large scale population studies, and only do so here to remoye 
what appears as somewhat incompatible Tea measurements of 
a few stars in our sample. As mentioned in Section |2] oscil¬ 
lations of stars with t'max ^10 /iHz cannot be reliably char¬ 
acterized with 80-day time series - the typical length of a K2 
camp aign. Prom the ^'max ex; relation (iBrown et alJ 

119911) . this essentially translates into a lower limit on log g, 
as indicated by the upper horizontal dashed line (Pig. [^). 
Similarly, the cadence of the data (^ 29.4 min) results in a 
Nyquist frequency of about 283/rHz, which defines an upper 
limit on t'max of about 270/iHz, and hence on logp (Pig.[^, 
lower horizontal dashed line). Oscillation frequencies aboye 
this limit will be too close to the Nyquist frequency, com- 
pro mising automated robust measurement of both r'max and 
Av (IStello et alJ2013l) . This includes stars oscillating beyond 
the Nyquist frequency. Again, our K2 results confirm these 
boundaries, with all but three detected oscillating stars falling 
within the APOGEE log g range 2.1-3.35. 

Comparing the APOGEE log g with the asteroseismic yal- 
ues, we see in Fig. [Jr that there is an offset of about 0.2 
dex between the two , in agreement with the findings of 
iHoltzman et al.l (l2015h . Note that we deriyed the seismic 
logp using the aboye iz^ax scaling relation with from 
APOGEE. We also note a large scatter, and in some cases a 
yery large deyiation (stars indicated by arrows). The seismic 
logp has a typical internal uncertainty of abo ut 0.03 dex for 
the length of data used here (IHuber et al.l2010l) . much smaller 
than can be obtained from spectroscopy. Hence we attribute 
the scatter in Pig. [3}) to the uncertainties in the spectroscopic 
determinations. This suggests an RMS s catter of 0.2 dex of 
the spe ctroscopic logp, in agreement with IPinsonneault et alJ 
(I20l4 : for a few stars a deyiation of up to 1.5 dex is seen. 
The latter extreme cases could potentially be blends in the K2 
data, where a more eyolyed, intrinsically brighter star with 
lower log g is detected seismically, while a less eyolyed star is 
the source in the APOGEE spectra. If they are not blends, the 
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Fi g. 3. — (a) los q versus Tefj from APOGEE IMaiewski et alJ2010l) for all 117 obsei'ved stars in our sample. We plot the values called ‘raw’ in Equations 2 and 
3 of lHoltzman et alH201.'5ll because their ‘corr’(ected) values are not available for all stars in our sample. Filled blue circles mark stars with detected oscillations, 
filled cyan circles show marginal detections, and empty circles are stars with no detected oscillations. A typical error bar is shown at (log g,Teff)=(3. 0,4100 K). 
Dashed lines mark the log g range of stars with izmax between 10/rHz and 270/rHz - the detectable range using K2 long cadence data s panning one full campaig n 
80 days). The vertical dotted line marks the upper limit in Tefj typically found for oscillating red giants in previous Kepler data iPinsonneault et ai]l2014ll . 
The gray solid curve shows a representative MESA stellar evolution track of a 1.2Mq model. Arrows indicate stars that potentially have a log g very different 
from the APOGEE value, (b) Relation between APOGEE and seismic logg of the stars with clear seismic detection. The solid line shows a 1-to-l relationship. 


two hottest stars marked by arrows must also be much cooler 
given their seismic signal, which indicates they are evolved 
red giants more luminous than the red clump. With these 
uncertainties in mind, we would expect that logp and Teff 
in Fig. for some stars are not necessarily representative 
of their true values. Hence, even in the scenario where we 
have 100% detection rates, we should find some stars with 
detected oscillations outside the log g-TgS region of expected 
detectable oscillations (Fig. [^, middle-left region bracketed 
by the dashed and dotted lines), and some non-detections in¬ 
side this region; this indeed appears to be the case. 

Finally, we examine our success rate in detecting oscilla¬ 
tions which, based on the above considerations, needs to be 
evaluated only for the stars expected to show oscillations. 
Within the log g-Tgg region of expected detectable oscilla¬ 
tions, there are 67 stars, of which 55 are clear detections, 7 
are marginal detections, and 5 are non-detections. We note 
that it is not inconceivable that all 5 non-detections have true 
values of log g and Tes that would make them fall outside 
the detectable region given that we find a similar number 
of detections outside the detectable region. However, some 
of these stars could be genuinely non-o scillating red giants 
suppressed by stron g binary interactions (iDerekas et aLlBOlTt 
iGaulme et alJl201^ . 

To judge whether the non-detections are simply caused by 
poor photometry, we show the white noise level versus ap¬ 
parent brightness for all 67 stars in Fig. |4] We see no clear 
trend between the detection category and the apparent mag¬ 
nitude, and hence noise level, which suggests that the non- 



Fig. 4.— Noise level in power spectra as function of Kp magnitude for 
all stars with APOGEE Tgff < 5000 K and APOGEE log g in the range 
2.1-3.35 (see Fig.[3i). The noise is measured as the median power in the 
range 260-280/iHz. For reference we show a fiducial line described by 
log(WhiteNoise) = 9.77 — 1.57-fCp + 0.0747fp^ ppm^//iHz that follow 
the lower envelope of the magnitude-dependent noise floor (red curve). The 
saturation limit of K2 is indicated by the dotted line. The increased noise for 
saturated stars can be mitigated using l arger aperture masks than in the VJ 
photometry used here ILund et alJ2Q15) . 
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detections are not generally caused by noisy photometry. We 
do see a possible hint of a faint limit at > 15 mag, but the 
low number of faint stars prevents any definitive conclusion. 
However, for the stars with detected oscillations, we calcu¬ 
lated the correlation between Kp and seismic log g and found 
it to be essentially zero {pKp,\o^g = 0.04). Hence, even our 
least-evolved giants (highest log g), which have the smallest 
oscillation amplitudes, span the entire magnitude range down 
to Kp = 14.75. Taking the uncertainty in Kp into account and 
being conservative, this result suggests that we will be able to 
detect oscillations in any red giant down to Kp ^ 14.5, as 
long as they are within the required log g range, which for the 
seismic logp is ~ 1.9 - 3.2. 

4. SUMMARY AND OUTLOOK 

We have performed initial asteroseismic analyzes of K2 
Cl data for over one hundred stars expected to be red gi¬ 
ants based on their spectroscopically determined \ogg and 
Teff. We detect oscillations in almost all stars cooler than 
5000 K v yithin 2.1-3.35 in log q (on the scale of the ‘raw’ 
data in iHoltzman et al.l l2015h . which comprise the target 
stars of ongoing K2-based galactic archaeology studies. The 
results indicate that our detection rate is complete down to 
Kp ~ 14.5, represented by the faintest oscillating star in our 
sample (Kp = 14.75) - a low-luminosity and hence low- 
intrinsic-amplitude giant with a detected i/max of 150/rHz. 
More stars fainter than Kp = 14.5 need to be analyzed in 
order to determine the true faint limit, which will set the limit 
for a magnitude-complete sample of red giants within the 1.9- 
3.2 range in seismic logp. Knowing this selection function is 
crucial for galactic archaeology. 

We note that the characteristic timescale for the spacecraft 
attitude control system to react to pointing errors was adjusted 
from C3 onwards to obtain lower pointing jitter (by about a 
factor of 3-4), which is expected to result in a total noise level 
only 30% higher than in the original Kepler mission (Doug 
Caldwell, private communication). We anticipate that this 
would push the faint end of the detection limit to even fainter 
stars than the Cl results reported here. 

Looking ahead, the prospects for characterizing oscillations 
in red giants using K2 data are very promising. We can ex¬ 
pect to be magnitude complete in our detection rates down to 
Kp ~ 14.5, and possibly fainter in future campaigns, which 
is strong affirmation for using red giants as distant probes of 
the Galaxy’s structure and evolution out to at least ^ 7kpc. 
Of order 5000-10000 targets have already been observed dur¬ 
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